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Abstract 

Recent developments in string theory have brought forth a considerable interest in 
time-dependent hair on extended objects. This novel new hair is typically character- 
ized by a wave profile along the horizon and angular momentum quantum numbers 
/, m in the transverse space. In this work, we present an extensive treatment of such 
oscillating black objects, focusing on their geometric properties. We first give a the- 
orem of purely geometric nature, stating that such wavy hair cannot be detected by 
any scalar invariant built out of the curvature and/or matter fields. However, we show 
that the tidal forces detected by an infalling observer diverge at the 'horizon' of a 
black string superposed with a vibration in any mode with / > 1. The same argument 
applied to longitudinal {I = 0) waves detects only finite tidal forces. We also provide 
an example with a manifestly smooth metric, proving that at least a certain class of 
these longitudinal waves have regular horizons. 
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1 Introduction 



One of the most confounding puzzles about black holes in General Relativity is the apparent 
incompatibility between their extremely simple structure, governed by the famous no-hair 
uniqueness theorems [|l|, and their generically large thermodynamic entropy 0. The former 
is indicating that the complexity associated with the latter is not encoded in the classical 
solutions. The crux of the no-hair theorems in General Relativity is that regardless of 
how a black hole was formed it is completely described once its mass, spin and Abelian 
electric/magnetic charges are determined. All other types of long range interactions, such 
as scalars with or without self-interactions, massive or non- Abelian gauge fields, etcetera are 
excluded because they would either lead to naked singularities ^ or would be unstable and 
thus would disappear shortly after the formation of the black hole More specifically, 
the machinery of the no-hair theorems rests on the assumption of regularity of a stationary 
null surface, i.e., the event horizon, and the demonstration that it leads to the vanishing of 
all matter charges, other than Abelian, which could carry long range interactions. Because 
the Equivalence Principle demands that the fields be coupled minimally to gravity, the 
vanishing of charges in stationary spherically symmetric space-times implies the vanishing 
of the external fields or hair. A small deviation from these results are black holes with non- 
Abelian hair, where the hair arises because of the nonlinearities in the gauge field equations 
of motion; however, this hair is unstable to perturbations, with the hairy black hole rapidly 
decaying into a bald one 0. The final result is that all of the structure of the configuration 
which collapsed into a black hole not encompassed by mass, spin and Abelian charges must 
have been radiated away in the process of formation of the hole. 

Recent developments in non-minimal models of gravity, and especially in string theory, 
where the modifications of General Relativity have a firmer foundation, have brought about 
some explicit counterexamples to the no-hair theorems |Q. Namely, the models considered 
involve new non-minimal couplings which provide extra sources in the equations of motion. 
However, whereas the particular solutions to these modified equations carry non-trivial long- 
range hair exterior to black holes, there are no new constants of integrations and so these 
aberrations of the no-hair theorems are deviations in letter only and not in spirit. The new 
charges are completely determined by the charges already present in Einstein's theory, and 
the non-minimal hair is secondary (as opposed to the primary hair carried directly by the 
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black hole charges) .0 Nevertheless, the secondary hair does affect thermodynamic properties 
of black holes, changing the expressions for both the temperature and the entropy [Q. 

Another kind of black hole hair may appear if we abandon the constraint of stationarity 
0-|Tl|]. Recent developments in string theory have brought forth a considerable interest in 
time-dependent hair on black holes and extended objects. In particular it has been suggested 



that such hair might give a classical accounting of the black hole entropy |TT|, |T5|. In any 
event, one might expect that this situation would be a more realistic description of a black 
hole, because stationary black holes could exist only as hermits in complete solitude; yet their 
very nature precludes this, as their own gravitational field permeates the whole Universe and 
communicates with all of its inhabitants. The principal obstacle to studying non-stationary 
problems is an enormous complexity of gravitational equations of motion, often resulting 
in their intractability. Indeed, there are comparatively few exact non-stationary localized 
solutions known to date, and none of them corresponding in all important aspects to the 
physical picture of black holes. Still, some of these solutions may provide useful models 
to study realistic hairy black holes. Most of these solutions arise because we can embed 
four-dimensional black holes in higher dimensions, via a procedure inverse to dimensional 



reduction (thus, the term oxidation — see e.g., ■'■^ certain special cases, the resulting 

oxidized solution has a null, hypersurface-orthogonal. Killing vector [|l3^[|H]- They can be 

used as a natural starting point for the construction of a more general family of non-stationary 

solutions characterized by a set of arbitrary functions. The presence of this isometry indicates 

that there exists a much larger family of wave-like solutions which reduces to the black 

hole-like solution for the special choice of its degrees of freedom. A similar correspondence 

exists, for example, between the Brinkmann wave solutions in General Relativity and the fiat 

space. We can get a glimpse at this larger family by using a solution-generating technique 

defined by Garfinkle and Vachaspati as a method to restore some of the original wave 

degrees of freedom. The resulting solutions still possess a null hypersurface orthogonal Killing 

vector, and thus are describing a gravitational disturbance propagating through the original 

environment at a speed of light i.e., a gravitational wave. 

A fundamental question of interest is to determine if this larger family still has a regular 

null surface, which can be reached by causal or null geodesies, i.e., an event horizon. In this 

^The secondary hair will be stable because the monopole charges do not vanish, in contrast to the non- 
Abelian hair mentioned above. 
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paper, we will examine this question for a certain family of five-dimensional black strings 
with a single wave profile function in an (/, m) multipole mode. First, however, we will give a 
theorem of purely geometric nature (valid in any theory of gravity which assumes a Pseudo- 
Riemannian geometry as the basis for the description of gravitational interactions) that such 
wavy hair generated by the Garfinkle-Vachaspati technique cannot be detected by any 
scalar invariant built out of the curvature and/or matter fields. However in our example, we 
will show that the tidal forces measured by an infalling observer diverge at the 'horizon' of 
the black string superposed with a vibration in any mode with / > 1. Hence the solutions 
with excited dipole or higher multipole modes contain null singularities. The same argument 
applied to the monopole mode shows that the leading order tidal forces are finite. We will 
also construct a class of monopole wave profiles for which the metric is manifestly smooth. 
Hence at least some of the wavy strings indeed have regular horizons. 

The paper is organized as follows. In the next section, we give a review of the wave 
generating technique, outlining the conditions the matter distribution must satisfy for the 
method to work. In section 3 we present our theorem on the elusive nature of Garfinkle- 
Vachaspati waves, and give its detailed proof. Section 4 contains the derivation of the explicit 
form of the wave-black string solution and a brief discussion of some of its properties. We 
show that these solutions contain null singularities for all higher multipole modes (/ > 1) later 
in section 4. In the last section, we construct an explicit family of longitudinal wave (Z = 0) 
solutions for which the metric which is analytic on the horizon (and not only continuous, 
as shown in 0). Finally, we present our conclusions and consider directions for future 
investigations. 



2 Garfinkle-Vachaspati Waves 

In this section, we review the wave-generating technique of Garfinkle and Vachaspati [|T7|. 
In doing so, we outline the situations in which it may be applied, i.e., the restrictions which 
must be satisfied by the metric and matter fields. The crucial requirement is that the 
solution possess a null, hypersurf ace-orthogonal. Killing vector. The resulting presence of a 
null coordinate allows one to effectively 'linearize' Einstein's equations, and restore hidden 
wave degrees of freedom. The method was originally proposed for the Yang-Mills-Higgs 
system coupled to gravity, and applied to straight cosmic string solutions in four dimensions 
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T7| . It was later extended to gravity coupled to a scalar and a two-form potential in five 



dimensions in the context of low energy string theory |T^. It was also used to construct 
similar wavy axionic strings in refs. [P, The following description will consider gravity in 
arbitrary dimensions coupled to a matter sector including various scalars as well as a set of 
different p-form potentials. Thus this discussion shows that this technique is quite generally 



applicable in supergravity or low-energy string theories |19 
Let us assume an action of the following form: 



/ = / d^'x^ I^R{g) --Y.Umct>af - 2 E/p(</')^{P+i)j • (1) 

Thus as well as the metric, we have included a collection of scalar fields (pa-, which appear 
with non-derivative couplings in the coefficient functions, /?.a(0) and fp{(f))- The above action 
also involves a set of p-form potentials, ^(p), through their field strengths: = dA(^py 

Hence there is a(n Abelian) gauge invariance associated with these fields, 5A(p) = (iA(p_i). 
Note that this action is written in terms of the Einstein-frame metric, by which we mean 
that there are no couplings to any (pa appearing in the Ricci scalar term of eq. (||). Now the 
gravity equations of motion may be written as 

i?^ - ^6^,R = ^T^ (2) 

where 

while the matter field equations may be written as 

= d,{^Ki^)g^''dM - lV9T.^W^^W,<Pa (4) 

b ^Va 

= d,(^jUmFi,^i)T""") ■ (5) 

For later purposes, we have been explicit about the appearance of the spacetime indices in 
these equations. 
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Now within this theory, let us assume a solution for which there exists a 



vector field k'^ which is: 

null : PA;^ = , (6) 

hypersurface orthogonal : V[^ku] = k[^'Vu]S , (7) 

and Killing : V(^A;j.) = . (8) 



Combining these equations, it is easy to show that k has a vanishing Lie-derivative on S, 
i.e., CkS = k^dfj_S = 0. Since we wish k to yield an invariance of the full solution, it is also 
assumed to have a vanishing Lie-derivative on the matter fields, i.e., the matter fields are 
form-invariant along the fiow of k. Hence, 

Ck(pa = fc'^9^0a = O (9) 

£fcF(p+i) = {dik + ikd)F(^p+i) = dikF^p+i) = (10) 

where the latter uses = di^ + ivd on forms with iy denoting the interior product, and the 
Bianchi identities, = 0. This form- invariance of the matter fields guarantees that 

the stress-energy tensor is form- invariant, as it must be given that Einstein's equations (0) 
are satisfied. Note that for the form fields, the vanishing Lie derivative is imposed on the 
physical field strengths rather than the gauge-variant potentials — the latter may vary by 
a gauge transformation along the A;- fiow, £a,.A(p) = (iA(p_i). In the following, we will further 
require that these fields satisfy an additional transversality constraint, namely, 

ikF(p+i) = kA 6'(p_i) . (11) 

Here, the right-hand-side is the wedge product of one-form k^dx^ with some (p-l)-form 
0(p-i), which necessarily satisfies ikO(p-i) = since (ifc)^F(p+i) = — see section 3.1, for 
details. 

The solutions satisfying the above conditions can be interpreted as gravity waves. Con- 
sider the coordinate system adapted to the fiow of k — as well as the cyclic coordinate v, 
there is a coordinate u, given 'roughly' by the integral of the dual one-form k = k^dx^. The 
vanishing Lie derivatives, eqs. (||), (|) and ([T0|) , simply means that none of the fields depend 
on V. Hence the only 'time' dependence can arise through the null coordinate m, and hence 
represents perturbations moving at the speed of light in a certain direction of the space- 
time. The Garfinkle-Vachaspati (GV) solution-generating technique extends the solution by 
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restoring additional wave degrees of freedom. To do this, they define a new metric by [|T7 



g'^, = Qf.u + e^'^k^K, (12) 

while leaving all of the matter fields unchanged^. The configuration {g', (j)a, Ap) also yields 
a solution provided the function \E' satisfies appropriate constraints. These restrictions will 
guarantee that after the metric is modified the form of the explicit fields appearing in the 
equations of motion remains unchanged. Hence eq. (12) is promoted to a map on the space 
of solutions. In order to maintain the wave interpretation of the new solution, first one 
requires that k has a vanishing Lie-derivative on \E', i.e., k^d^^ = 0. One may verify that 
this ensures that the hypersurface orthogonal and Killing conditions, eqs. (|^) and (^), are 
still satisfied with the new metric (with the same S). It is also obvious that the null condition 
(^) still holds with g'. To determine what other restrictions must be imposed on \I', one must 
consider the changes which the map ([T^) induces in the equations of motion, eqs. d^-^l). 



We begin by demonstrating that the matter field equations are invariant under ([1^) . 
First note that the determinant of the metric would only be modified by terms proportional 
to k^^kfj^. However since k is null, the determinant is invariant: det{g') = det{g). On the 
other hand, the inverse to the metric g'^^^ is given by 

g'"" = g^"" ~ e'^ <i/ k^k" . (13) 

Now given the transversality constraint (^Ij) , we see that raising all of the indices of the form 
fields with the new metric yields 

where in the first line, terms with more powers of k have automatically vanished since 

ikO(j,-i) = 0. Given this result as well as the invariance of det{g), it is clear that the 

equations of motion for the form fields (^ are left unchanged by the map (|T^). Further 

given the vanishing Lie derivative (^, it is clear that g'^^d^ipa = g^"dv<Pa and hence the scalar 

equations (H) are also unchanged. Thus we have shown that the configuration {g' , (pa, Ap) 

provides a solution of all of the matter field equations. 

^Throughout this section, A:^ = g^j/fc''. However, from the definition of the metric g'^^, and the fact that 
k is null, it doesn't matter which metric we use to lower and raise the index of k. 
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The same considerations as above also show that the stress-energy tensor with mixed 
indices T^^^, given in eq. @, remains unchanged when the metric is modified as in eq. ([T^) . 
The only step which remains in examining the gravity equation is to compute the change 
of the mixed-index Einstein tensor. This is reduced to computing the change in the Ricci 
tensor only, because the Ricci scalar may be eliminated from eq. using R = -^^m ^ 
dimensions. It is straightforward to calculate that the Christofel symbols for the two metrics 
are related by T'^^^ = T'^^ + Q'^^^ where 

^'.x = ^(V,(e^vl>Pfc;,) + VA(e^^FA;,) - V^(e^vl>r A;;,)) . (15) 

Further using the properties of k, as well as k'^d^S = = k^d^J^ , one can show that 
= = kt'Vx^^u = ^U^pu- Thus one finds that R'x^ = Rxu + ^ p^iy and hence 

B!\ = i?^ + e^^fc'^fc^i^A. + Vp(fit(7'^^) (16) 

Again using the properties of /c, one can show: k^Rxu = ky\/'^S/2. Substituting this into 
eq. (p!6D and after a few more manipulations, we finally arrive at 



R'^'^ = _ le^Pfc^V^^ (17) 

Therefore, the variation of the mixed Ricci tensor R^^ under (0) is proportional to V^\l', and 
so vanishes if we demand that \1/ solves the covariant Laplace equation (in the background 
defined by g^^,). 

Therefore given a solution for which eqs. (p|-pT]) are satisfied, then eq. (|T2p provides a 
map to a new solution provided 

= and V^^ = . (18) 

We note that the condition V^\l/ = is really just the spatial Laplace equation since k^d^'^ = 
0. Furthermore the latter indicates that the moduli of the new solutions only depend the 
retarded time - u, above - and therefore they still represent gravitational waves. 

The preceding discussion was phrased in terms of the Einstein-frame metric. However, 
in some instances (such as the example in section it is more convenient to work in terms 
of a conformally related metric, e.g., 

gpu = e-^^^gp. . (19) 
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Given in terms of g, the action will contain non-minimal couplings of the scalar fields to the 
Ricci scalar. It is straightforward to show that most of the constraints (l^^lTD are unchanged 
when written in terms of the conformally transformed metric. The only change is that the 
hypersurface orthogonal condition becomes 

^[fiK] = k[pyu]S where S = S — a(0) . (20) 

Here we have denoted = g^^k^ . Thus the map (|T2|) becomes 

~g'^,, = ~g^.u + Q^^~k^,k. (21) 

which is equivalent to g'^^ = e^^'^^g'^^- Finally, of course, the constraints on which ensure 
that eq. (0) provides a solution of the gravitational equations of motion are identical to 
those appearing in eq. ([18|) . The latter may be written as follows, when expressed in terms 
of g: 

k^'df,-^ = and df,(^e^''''^^^f^g'"'d^^^ =0 (22) 

where again D denotes the dimension of the spacetime. 

Before concluding this section, let us mention a few simple extensions of the original 
action (|l]) for which the preceding discussion would still be applicable. First, we could add 
a scalar potential U{(f>) (which could then include a cosmological constant) to eq. (|l]). This 
would modify the stress-energy (^ by a term proportional to 6^,yU{(f)), as well as adding a 
d^^U term to the scalar equations (^. Both of these terms are obviously invariant under 
the map (^2|) and so the above construction is still valid. 

Another non-minimal coupling, which commonly arises among the form fields, are Chern- 
Simons-like terms appearing in the definition of the field strengths, e.g., for some choice of 
m and n, F(^rn+n+i) = dA{m+n) + tt^{m) A dA(^n) where we still assume F(^m+i) = dAf^^n) and 
-F(„_|_i) = dAi^n)- This definition leads to two modifications in the equations of motion. (Note 
that for F(m+n+i) the Bianchi identity is also modified, dFi^m+n+i) = ct-F'cm+i) A -F(n+i) , but 
this has no consequences for the present discussion.) First, eq. (^) is modified for p = m by 
the introduction of a source term which is proportional to 

With the same transversality constraint on F(m+n+i) as above, eq. ([TTl) , we still conclude that 
its form with all indices raised remains invariant, as in eq. ([I^), and hence this new term is 



8 



also unchanged by eq. (|^) . The second change is the appearance of a source term in eq. (|) 
for p = n. In this case upon applying the equation of motion (|^) for p = n + m, this second 
source term takes precisely the same form as above, merely interchanging the roles of m and 
n. Therefore it is also invariant under the map (p!2D. Thus the GV procedure will still be 
valid for solutions where these Chern-Simons-like couplings make nontrivial contributions. 
The final extension which we will consider is the addition of topological interactions, e.g., 

where i + m + n + 2 = D. Such an interaction again modifies the equations of motion (§) 
for p = i,m,n hj the introduction of source terms of the form, e.g., 

where e is the Levi-Civita tensor in D dimensions. Given the invariance of the determinant 
of the metric, e' = e as forms. Hence we must compare raising all of the indices with g' and 
g. Given eq. (0), we see 

= e^^-P (24) 

where the vanishing of the second term in the first line relies on a result proven in the 
following section — see eq. (|40|) . Therefore the new source terms are again unchanged by 
the map (|12D , and so these topological interactions provide no obstruction for this solution 
generating technique. 

3 The Elusive Wave 

Generically it is the case that the original solution and that carrying a wave induced by 
eq. (|1^) are not diffeomorphic, however, the wave turns out to be very elusive. In this 
section, we will present a theorem of a very general nature, showing that all the scalar 
curvature invariants of the two metrics related by the wave-generating technique are in fact 
identical. Thus, no curvature invariant of the oscillating metric can be used to detect the 
presence of the gravitational wave! The theorem generalizes in a straightforward way to 
include any scalar invariants constructed from both the metric and matter fields. 
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To prove this result, we need only assume the existence of a metric g^,y which admits 
a null, hypersurface-orthogonal. Killing vector. We do not require that the metric solves 
Einstein's or any other equations. Thus our result is purely geometrical in nature, and holds 
for any metric satisfying the symmetry condition. The precise statement of the theorem is 
as follows: 

If g^u is a pseudo-Riemannian metric admitting a null, hypersurface-orthogonal, Killing 
vector , and g'^^ = g^u + i^k^k^ where k, is any scalar Lie-derived by k to zero, i.e., C^k, = 0, 
then all of the scalar curvature invariants of g'^^ are exactly identical to the corresponding 
curvature invariants of g^^. 

Hence the GV transform ( ]T2| ) provides one example of such metrics with k = e'^\&. In 
the following, we will refer to g^^, as the original metric, g'^^^ as simply the primed or shifted 
metric, and k as the wave profile, in analogy to the last section. 

Before presenting the details, let us first give a brief sketch of our proof. Before examining 
the curvature invariants, we establish two crucial results. First, contracting k^ with any 
tensor constructed from the original metric, its curvature and covariant derivatives of the 
curvature and/or any scalars with a vanishing Lie-derivative under k, produces a sum of 
terms in each of which k appears uncontracted. Second, any tensor {e.g., the curvature or 
covariant derivatives thereof) in the primed background may be written as the sum of that 
for the original metric plus a /t-dependent term, for which all of the contributions are at 
least bilinear in the Killing vector k. These two results set the stage for an examination of 
the scalar curvature invariants. There we find that all of the new ^-dependent terms vanish 
by combining the previous two results with the fact that k is null. Hence, we conclude 
that the original and the corresponding primed invariants are identical. As a consequence, 
no evidence of the wave profile k can be detected in any of the curvature invariants of the 
metric. This result represents a generalization of the previous studies of curvature invariants 
of geometries admitting null, covariantly constant. Killing vectors [^, An immediate 
corollary of this theorem is that if g^^ represents an extended black object with a regular 
horizon to which we add a GV wave as in the last section, then the oscillations will not 
produce new scalar curvature singularities — however, in later sections, we will discuss the 
limitations of this statement with explicit examples. 

0) Useful Formulae 

Let us start by listing some important formulae which arise from the existence of a 
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null, hypersurface-orthogonal, Killing vector. As stated before (in eq. (J^)), the hypersurface- 
orthogonal condition amounts to 'W[^ky^ = k[^'Vu]S for some scalar S which can be determined 
from the metric. Combining this with the Killing condition (§) yields 



V^K = -{k^V,S-KV^S) . (25) 



Further, it is not difficult to see that the Killing condition alone leads to 

V.V^A;a = kPRp^^x (26) 

as can be determined by considering the commutator of two covariant derivatives acting on 
k. We also note that we may express 

= k[xV„]V,S~l-k[xV^]SV,S (27) 



Ul...Uq- 



using eq. (pS]). Recall the expression for the Lie-derivative of a general tensor T/^i - Aip 

^V-'- U\...Vq '-' * A-^ Ul...Uq 

-T^-^\,...uqV xv^' - ... - T^^-\i...., Vat;'^'' (28) 
+T^''-^- x...uqy + ■■■ + T^''-^\^,„xS/ ,y . 

From this definition and the Killing condition (^, one may show that the Lie-derivative with 
respect to a Killing vector Ck commutes with the covariant derivative. Begin by considering: 

^^j0.kTu,...uq = k^V^VxT,,...,^ + VaT,,....,V^A;^ + . . . + V^T,„„xV,^k^ 
+Tx...uqVf,V,y + ... + V^T,,...AV^V.,fc^ 
= A;^VaV^T,,...., + VaT,,...,, V^A;^ + . . . + VpT,„„xV,^k^ 

~\~(Tx..Mq{R Vipfi "I" RpfJ.Ul )k^ -|- . . . + T,^j^ x[R UqPfl + RpilVq ) 

= CkV^T,,,,,,^ (29) 



using eq. (Pq) and the standard commutator: ^ x\T^^y^ = R^'uiXpTp. + . . . . Hence, 
V^]T^i...i,g = 0. The case where V^] acts on a tensor with some raised indices is 
trivially related to this one because Ckg'^^ = = V ^g^^. Similarly, one can show CkRpXvp = 
0, as well as CkS = /c^V^S* = 0. Further recall by definition, Cki^ = k'^V p,K, = 0. We will 
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find all of these formulae useful below, when evaluating the contractions of tensors with the 
null vector k. 

1) Contraction Identities 

The first step is to show that for any tensor built out of the original curvature, any 
scalars with vanishing Lie- derivative under k {i.e., S and k), and any number of covariant 
derivatives (with respect to the original connection) acting on either of these, the contraction 
k^T^i...upfiXi...Xq is at least linear in vector k^, i.e., it can be expressed as a linear combination 
of terms which factorize as some tensor of rank lower by two and the vector k: 



k^T^l..MpfJ.\l...\q — X! k'^n^ii...Un..Mp\l...\g + X! 6^^^ . . .jy^ Ai . . . A^. . . A, (^0) 

n=l n=l 



where underlining an index denotes its deletion from the expression. (For convenience, we 
will work with only covariant, i.e., 'downstairs', indices.) This is true in the few simple cases 
encountered so far, e.g., k^V^S = = k'^V^K and kPRpyx„ as given in eq. (^). Let us now 
show that it holds in general for what we denote as 'primary' tensors, namely tensors obtained 
by an arbitrary number of covariant derivatives acting on the curvature or the scalars S", k. 
The proof again relies on mathematical induction, and makes essential use of eq. (^) which 
allows any covariant derivative of k to be re-expressed in terms on undifferentiated fc's. First, 
we establish the result for the simplest cases: for any scalar B Lie-derived to zero by fc, we 
have k'^'V/^iB = and further 



1 
2 



--k^V^'SV^B (31) 



using eq. (pS]). Now combining CkRap-ta = and eq. (pTf), we find 

k^V fj.Raf3^cT = R'^ f3-yak[fjy a]S + . . . + Ra/S-y'^ k^^V a]S 

qW , , L -9(4) 



kJZ + ■■■ + kj± (32) 



where in the second line we collect the like terms. Explicitly, one finds, e.g., 6^^^^ = 
-\{V^^R^p^. + V/3V[^5 V<,]5) and ^(2), = -\{V^R^^^. - V„V[-,5 V.]5). Similarly, 

k ^ ^Ral3"/cr = V^(A; Rap-ya) ~ Ra/S-ycr^ f^k 

= k,e^^l + ... + kj^^l (33) 
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again using eqs. (|^) and (pTl) . 

So now let us assume that eq. (|30|) holds for all primary tensors of rank q or less. Now by 
our definition, a primary tensor of rank g + 1 will be obtained by covariant derivative acting 
on a primary tensor of rank g, i.e., V \Tx^ Xq- Hence we consider 

= E Va, {kxAl>^,„x) + ^T,a,...a,PVa,5 - \T,x2...Xqkxy'S 

n=2 ^ ^ 

= E (A:a„ Va,<La,...A, + VA,fcA„C^A,...A, + ^/^A„<^..A,...A Va,^) 
n=2 ^ Z / 

-^A^A^T^A^.-AgV^^ (34) 



E ^A„6'aiLa„...A„ 



n=l 



after collecting the like terms into the tensors 0^'^\ Of course, the same result follows if k^ 
is contracted with any of the other indices on T above. The last possibility is k^V pT\^,,,\^. 
Here we use CkT\^,,,\^ = which holds by = and the fact that both S and the 

curvature have a vanishing Lie-derivative under k. Thus, 

fc''V^TA„..A, = -T^A2...A,VA,fc'^-...-TA,..A,_,^VA,fc'' 

= ^(t^A2...a,Va,S + ... + Ta,..a,_,^Va,5)A;'^ (35) 



^ ( Wa,A;a, + . . . + Ta,...a,.,^A;aJ 
E kxM^. 

n=l 



where we applied eq. ( pil| ) for rank q primary tensors in the second line and gathered the 
like terms. With eqs. (|3^ and (|35|) , we have shown that eq. ( ^0|) holds for any index on a 
primary tensor of rank g + 1. Therefore by induction this factorization property is established 
for all primary tensors. 

At this stage, we must consider the 'secondary' tensors, i.e., the 9 tensors produced in 
the contractions of k with primary tensors — although here we will leave many of the details 
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to the reader. Considering the simplest examples in eqs. ( pTf ) and (pI H33|) , one finds that 
the 9 tensors are not simply primary tensors, but rather involve certain products and/or 
contractions of primary tensors. However, in those particular examples, it is not hard to 
show that they share two important properties in common with the primary tensors: i) the 
^'s have a vanishing Lie-derivative under and ii) upon contraction with fc, they factorize 
as in eq. ([30|). Having established that these conditions apply in the simplest cases, it 
is straightforward to formulate an inductive proof to show that they also apply for the 9 
tensors produced from primary tensors of higher rank. One would begin by assuming that 
(i) and (ii) hold for the ^'s arising from rank q primary tensors, and then examine those 
produced at rank g + 1. In the case considered in eq. (|33), one has 



_ \ 1 \ 

^ n=l ~ ^ 

Ca,...a, = V.,^£ forn>l. (36) 

It is easy to show that 9^^^ will satisfy both the Lie-derivative and factorization properties 
because all of the components of which it is comprised {i.e., 9^"'^^s, T, VS") do. For n > 1, 
^(") will be Lie-derived to zero by k because ^("^'s are and [Ck, V^] = 0. Given the vanishing 
Lie-derivative and that factorization (130|) holds for the ^^'■"^'s, one would extend this condition 
to the ^'•"^ tensors in the same way as in eqs. (|3^) and (|35|) for the primary tensors. One 
should then examine the ^'s arising in eq. ( pS]) for that particular case of the rank {q + 1) 
primary tensors. One is again able to show that both (i) and (ii) applies for these 9 tensors 
as well, although the index gymnastics is somewhat more involved. Furthermore, the proof 
of these properties extends in a similar way to any 'higher-order' tensors, that is any new 
9^s produced by contracting k with 9 tensors. 

Now this intermediate result for the 9 tensors is necessary in order to show that the 
factorization property (^) also holds for products of primary tensors, with arbitrary con- 
tractions of pairs of indices. Our principal tool here is again mathematical induction. We 
begin by considering quantities of the simple product T^ll,,yp^j^Xi...Xg'^'^'^''^^"'^''i^i---'-jp where both 
T(i) and T^^) are primary tensors. Using the properties of primary tensors, we find 

p 

r,/^7n(l) ^(2)Ai...A, _ r n{n,l) y.(2)Ai...A, 

lyl...^l...UpXl...Xq'^ UJ1...UJI — ^ '^iyn"l/l...l^...l/pXl...Xq'^ K'^ ' J 

n=l 

q 

V A'^ rp{2)X^...Xq 
^ '^Xn V\...VpX\...\n...Xq^ 1^1. ..W; 

n=l 
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(38) 



Here the first term has the required form, and so for q = 0, i.e., no contractions, we would 
have the desired result. However, for g > 1, we need to look at the cross-terms in the second 
sum e.g., 

q 

,2)X2...Xn...Xq 

'^^l^Ul...UpX2...Xg-^ COl-.-Ul — Oui...Vj,X2...Xq 2^'^ ^ U11...UJI 



n=2 

,n(P+hl) u n{q+n,2)Xi...Xg /oq\ 

n=l 



Here we will have arrived at the desired form if g = 1, but for g > 2, we have generated further 
cross-terms in which k is contracted back on However, since the ^'s also factorize 

according to eq. (^) as described above, we may continue this procedure. Now given that q, 
the number of contractions, is finite, this 'ladder' of k contractions will eventually terminate, 
since at each step the number of contracted index pairs is reduced by one in each of the 
subsequent cross terms. Therefore after a finite number of steps, we arrive at a factorization 
without any contractions yielding the desired form. Given this result, it is obvious that we 
may in a similar way extend eq. (|30D to apply for an arbitrary product of primary tensors, 
including arbitrary contractions. 

To conclude this subsection, we consider contractions of k with the Levi-Civita tensor, 
e, i.e., the volume form on the £)-dimensional spacetime. To analyze this case, we resort 
to local coordinate patches adapted to the properties of k. First, the Killing condition 
(§) indicates that we can find a cyclic coordinate such that k'^d^ = (d/dv). Next the 
hypersurface orthogonal condition (|^) indicates that we can find a dual coordinate such that 
k^dx'^ = e~^du. In this local coordinate patch, one of the free indices in k^e^a---i3 must then 
take the value m, and hence we can write 

ike = kAe (40) 

where in this case 6 is some (D-2)-form. However, as written this result is coordinate 
independent and so must hold in general.^ Thus we have shown that the Levi-Civita tensor 
factorizes in the same way as the primary tensors. 



Eq. ( [4 of ) is the essential result required to prove eq. (24) 
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In fact, the resulting 9 satisfies eq. (|30D as well in a trivial way since ikO = 0. This result 
is again most easily derived using the local coordinates introduced above. Because of the 
antisymmetry of the Levi-Civita tensor, none of the indices on the right-hand-side of eq. ( ^OD 
correspond to the cyclic coordinate v. Hence ikO = with this choice of coordinates, but 
this equation must then be valid in general. 

As an aside, we note that the preceding discussion is equally applicable for the ^(p-i) 
arising in the transversality constraint ([Tl|): = k A In this case using the 

adapted coordinates, the antisymmetry of the field strength -F{p+i) ensures that 6'(p_i) does 
not carry a v index. Hence one finds that ikO{p-i) = 0. 

Given the previous results, our final conclusion is that eq. ( pOD applies for an arbitrary 
product of primary tensors and Levi-Civita tensors, including arbitrary contractions. 

2) Shift in Tensors 

Next we consider the difference between tensors calculated for the original and primed 
metrics. First the shift in the connection coefficients Q'^x — ^'^x ~ ^ux given by 

^ux = l{v,{t^k^kx)+Vx{^k^k,)-V^{f^k,kx)) (41) 
= ^ {k^'kxV^K + kf'k^VxK - k^kxV^n) 

~ {k'^kxV.S + k'^k.VxS - 2k,kxW^S) (42) 



using eq. (P^l). The corresponding shift in the curvature is then R'^uXa = R^uXa + Va^^^o- ~" 
V a^^x + ^vo-^pA " ^ux^'pa- Howcvcr, the vanishing Lie derivatives, Cki^ = = CkS lead to 
^i^o-^pA ~ ^ux^'pa — 0- Hence the shifted curvature reduces to 

Bl^uXa = R^uXa + Va^^. " ^ a^uX (43) 

Hence the curvature with all covariant indices is R'^^xa — ^puXa + V^a^^Jj, + k^k^Rp^xa 
Now using eq. (|27|), as well as rewriting 2(7^pV[a^ct1i/ "with eq. (|25|) , we arrive at 



R'pu^a = Rf^uXa + 2A;[^V^] V[Afi: k^] - 4k A;[^V^] ViaS* k^] 

+3ki^Vu]nk[xV„]S + 3k[pyu]S /i;[aV^]K - duki^V^jS kixV„]S (44) 

Here the key observation is that the primed curvature can be decomposed as R'^^xa — 
RfifXa + XpuXu where the tensor Xp-uXu is bilinear in the vector A;^. 
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In fact, a similar decomposition holds for any covariant derivative of the curvature. 
To establish this result, first note that V^R'^^xa = Vp-R^^Aa + XpfiuXa where XpfiuXa = 
^pXpuXa - ^'^p^.KuXa " ^'^puR'pu.Xa " ^"^pX^uu^a " ^pa^uXu^- ^ is Straightforward to show 
that the tensor XpfiuXa is bilinear in the vector k. This is clear for Vx using eqs. ( p5D 
and (^4]), and for the flR terms using eq. (0). Now similarly the fix terms are quar- 
tic in k, but a closer examination shows that the sum of these terms vanishes. Then, 
by induction, we see that if Xpi...pnpyXa is at least bilinear in fc^, then Xppi...p„puXa must 
also be so. This is again straightforward by combining the above formulae in the def- 
inition of X- Xppi...pnpuXa = VpV'p^. . .V'p^R'^^x^ -VpVp,...Vp^R^uXa = VpXp,...pnpuXa + 

fi'^p^Xuj...p„puXa + • • • + ^ppi V^^ . . . Vp^RpuXa + ■ ■ ■■ Heucc, all of the new ^-dependent terms 
appearing in the primed curvature and covariant derivatives thereof, Xppi...p„puXcT, are at least 
bilinear in k. A detailed inspection shows that xpi...p„puXa contains a sum of terms of order 
k'^P with p = 1, 2, . . . , 1 + \n/2\, where \n/2\ denotes the integer part of n/2. The precise 
powers will not be important below, the key point being that they all begin at k"^. 

Finally we close this subsection by noting that the Levi-Civita tensor is invariant under 
the shift between the original and primed metrics, i.e., s'^.-.p = £p--p- This result follows 
since, as noted in section 0, the determinants of the two metrics are equal because k is null. 

3) Scalar Curvature Invariants 

We are now ready to consider the scalar curvature invariants which can be built for the 
shifted metric. The most general invariant will consist of an arbitrary product of curvatures, 
covariant derivatives thereof^, and Levi-Civita tensors, with their indices (all assumed to be 
covariant) contracted by the inverse metric g"^'^ . Hence a generic term is of the form 

N K M 

^' = UKi....i n4j...Ai,n^'"'' (45) 

j=l ^ k=l 1=1 

along with a rule for contracting the upper with the lower indices. Here Y^jLi Qj + DK = 2M 

where D is the dimension of the spacetime. In fact, one need only consider = or 1 since 

the product of two Levi-Civita tensors can be reduced to a sum of products of metric tensors. 

Now from the previous subsection, we know that all of the tensors in the first product 

can be decomposed as T^^ = T^^ -|- Xpi...fMn where Xpi...fMn is at least bilinear in the 

■^Here one might also include covariant derivatives of the scalar S which are implicitly geometric tensors 
derived from the original metric. Admitting these extra tensors would not change our conclusions. 
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vector k. Further, we have e' = e while the inverse metric is given by g'^'^ = g^^ — nk'^k'^. 
Hence we see that the invariant (145|) can be decomposed as follows: 

N , ^ K M 

^' = n + X........ n ^Aj...A^ 11(9'"'' - ^k'^^k^') (46) 

j=i ^ ^ k=i 1=1 

= i+j 

where X is the invariant of the same algebraic structure as I' but constructed for the original 
geometry. The difference then contains all of the information about the wave. Simply 
multiplying out the terms in (^61) , we may write 

J = ^A;^^r^...A;^«r»f,,.,...,,.^ (47) 

i=l 

where the tensors T are products of primary tensors and possibly Levi-Civita tensors, in- 
cluding contractions by g^". Given the results of subsection 2), we are assured that JT" is 
at least bilinear in k. Note that antisymmetry of indices, e.g., in e, may eliminate certain 
contributions above. However, from subsection 1), we know that the tensors T factorize as in 
eq. (|30D when contracted with k. Hence we conclude that J' oc fc^/c'^ = and so I' = I. Thus 
we see that any scalar curvature invariant is identical for the original and primed metrics, 
which concludes our proof of the theorem. 

We note that we can easily generalize our theorem to cover scalar invariants constructed 
using both the geometry and matter fields. As in the section ^, we consider a matter 
sector including various scalars (pa, and p-form potentials A(^py We also require that these 
fields satisfy the same constraints as there: they are form- invariant along the k flow, i.e., 
'C-k'Pa = = £fcF(p+i) as in eqs. @ and (|l^). The field strengths are transverse to the flow, 
i.e., = k A 6'(p_i) as in eq. (pT]). Recall we also have ikO(p-i) = 0. Given these 

results, one can further show that Ck0(p^i) = — this is easily shown by referring to the 
local coordinate patches introduced at the end of subsection 1). 

Let us then reconsider the contraction identities proved in subsection 1). Given the 
properties imposed on the matter fields above, it is straightforward to extend the discussion 
to include the scalars, the field strengths, and covariant derivatives of these, as primary 
tensors which satisfy the factorization equation (^0|). Similarly one may also show that all 
higher-order tensors and hence arbitrary products of the primary tensors satisfy the same 
factorization property. Next as in subsection 2), we consider the shift in tensors calculated 
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for the original and primed metrics. While the scalars and field strengths themselves are 
not affected by the shift in the metric, using eq. (^T]) as well as eq. (p5D , it is clear that the 
shifts in covariant derivatives of these fields will always be at least bilinear in the vector k. 
Hence both of the crucial results established for the curvature and its covariant derivatives 
are easily extended to the matter fields and their covariant derivatives. Thus the invariant 
(^) can be extended so that the T also include these latter fields, and the same final result 
still holds, i.e., the invariant is independent of the wave profile k. 

In summary then, we find the rather surprising result that all scalar invariants, involving 
any number of covariant derivatives of the curvature and/or matter fields, are identical for 
both the original and the shifted metrics. The essential requirement that the original metric 
has to satisfy is to have a null, hypersurface-orthogonal. Killing vector, which is supplemented 
with certain constraints on the matter fields as well. Thus no scalar invariant contains any 
information about the wave profile k, and hence to determine how the geometry has been 
modified, one must consider quantities such as tidal forces or non-local holonomies. We 
emphasize that our theorem is of a purely geometric nature, and holds for any theory of 
gravity in any number of dimensions, as long as it assumes a pseudo-Riemannian geometry 
as the basis of the description of gravitational phenomena. 



4 A Five-Dimensional Black String 

Our original motivation in this project was to investigate the properties of extended black 
objects in higher dimensions carrying time-dependent or wavy hair. The theorem of the 
previous section tells us that coordinate-invariant probes are inadequate to examine the 
properties of such undulating solutions constructed through the GV technique. Hence to 
study the smoothness of the horizon in the presence of a wave in the next section, we are 
lead to consider the existence of parallelly propagated curvature singularities. In order to 
do so, however, we must consider a concrete example. Thus in this section, we present a 
family of undulating black strings which are low energy solutions of heterotic string theory 
in five dimensions. We begin with a stationary solution with a null hypersurface-orthogonal 
Killing vector, which results from uplifting a four-dimensional solution first written down 



by Cvetic and Youm |14]. Then we apply the GV technique to generate oscillations on the 



string. Similar oscillations of singular strings were considered in refs. [0, ^, and of black 
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strings, in refs. p|, [15 



The low energy action for heterotic string theory in five dimensions includes the following 
terms 

1 = J (^R{G) + 4(V$)2 - - (Va)' - ^e^'^F^ - ^e'^-F^) (48) 

as well as the metric, we have included two scalars, the dilaton $ and a modulus field a; two 
gauge fields with F^^ = d^A^ — d^A^ and F^y = d^A,, — du-A^; and the Kalb-Ramond field 
with 

Hf,ux = d^B^x - -{A^F^x + A^Fux) + (cychc permutations). (49) 

The metric above is the so-called string-frame metric. The Einstein-frame metric would be 
given by 

g,, = e-^^/^G^, . (50) 

With the latter metric then, the dilaton coupling in the Einstein term is eliminated and the 
action becomes 



/ = /rf^Xv/^(i?(^)-^(V$)^-^e-«*/3i7^ 



However, we choose to present our solution in terms of the string-frame metric, which has a 
much simpler appearance in the present case. 

We will be interested in the following black string solution: The string frame metric is 

h h 

while the remaining fields are given by: 



ds^ = Ldu^ + ldudv + M (dr^ + r\de^ + sin^ 6 d(f)'^)) (52) 



B = — du A dv 
h 

A = -Pi cos 9d(p A = P2 cos 9 dcp 

e'" = e/k e^* = ki/h^ . (53) 



where we defined the following functions: 



k = 1 + y £ = l + y. (54) 
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Hence this configuration is specified by four different parameters. The solution could be 
simplified by setting all of these equal, however, we wish to illustrate that our results apply 
in the generic case. We will assume that all of the constants are positive, in order that our 
solution properly describe a black string with a horizon at r = 0. 
Considering the asymptotic metric, one has for large r 

ds"^ ^ du^ + 2 du dv + dr"^ + r'^dVL = -dt^ + dy'^ + dr'^ + r'^dVL (55) 

where y = u + v and t = v. Hence we should consider y as the spatial coordinate running 
parallel to the string, while t is the asymptotic time. Note that as r — > 0, grr — 
indicating the presence of a degenerate horizon. Near the horizon, the metric becomes 



ds'^ ~ -pz-du^ + — dudv + P1P2 
Q2 Q2 



^^]\dn^ 



(56) 



The solution has two Killing vectors which are of interest 

k''d^ = d,=dt + dy h''d^=du = dy . (57) 

(These are in addition to the standard rotational Killing vectors for 6 and 0.) The first of 
these is the null generator of the horizon. In this role, k has the rather unusual feature that 
it is null everywhere - not just at the horizon. Further, it is not given by dt in the asymptotic 
coordinates, rather we have k^d^ = = dt + dy. Here the dy contribution is related to the 
presence of linear motion along the y direction. The coefficient may be interpreted as the 
'horizon velocity', which in the present case is one, i.e., the speed of light. 
The null Killing vector k is hypersurface-orthogonal as well, satisfying 

V[f,k^] = k[^Vy] log h (58) 



where h is the same function defined in eq. (|5^ . Hence the metric admits the symmetry 
desired for the wave-generating technique. It is also straightforward to show that the matter 
fields satisfy the appropriate conditions (|^-[TT|). Hence following the discussion of section H 
we apply eq. (|21|) to define a new string-frame metric 

G'=G^, + h^k,k, (59) 
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where following eq. (|22|) , \& is chosen to satisfy 

k^Vf,"^ = 0, and df, (e-^'^V^G^'^d^'^) = . (60) 

The first condition is simply d^^^ = 0, i.e., \E' is independent of v. In the present case, the 
second condition reduces to V^p"^ = 0, i.e., Laplace's equation on a flat spatial metric in the 
transverse coordinates {r,6,(j)). Thus the general solution for eq. (^) may be written as 

* = E ^' + bimiu) r-('+')) YUO, 0) (61) 

l,m, 

where Yim{0, 0) are usual spherical harmonics, and aim and bim are arbitrary functions of u. 

Let us consider the various perturbations in turn. The case with and / = yields an 
asymptotically fiat metric, which in fact is diffeomorphic to the original, 

ds''^ = -(f + a(u)) du^ + -dudv + ki(dr^ + r^dQ) 
h h 

= -du"^ + -dudv + k£(dr'^ + r'^dn) (62) 
h h 

where v = v + ^a{u)du. Note that the constant term in / is a special constant case of these 
perturbations, and so could also be eliminated in the same way. With the choice r' and 
/ > 1, the metric is not asymptotically fiat, and so we do not consider these solutions as 
providing perturbations intrinsic to the black string, i.e., 'wavy' hair. Instead they would 
be more accurately described as embedding the string in a space filled with (asymptotic) 
gravitational radiation. The same is apparently true for the r' mode with / = 1, but in fact 
this solution yields an asymptotically fiat metric, as is seen as follows H, ||: First introduce 
Cartesian coordinates on the transverse space, in which case the wavy metric becomes 

1 2 

ds' = — (f + ai(u)x^) du^ + — dudv + kidxidx'^ . (63) 
h ^ ' h 

However the following coordinate transformation, 

= x'-A' (64) 

with 2Ai = Oj, produces a metric which is manifestly asymptotically fiat 

1 2 
dS'^ = T{f+ {hk£ - l)i2j du^ + -du (dv + {hki - l)Aid5:^ + kidxidx' . (65) 
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These waves represent oscillations of the string in the transverse space. 

The perturbations generated with r^(^+'^^ are all localized near the horizon, and leave 
the metric asymptotically fiat. Hence we may consider these deformations as candidates 
for 'wavy' hair on the black string. For / > 0, these deformations produce g^u — > C)0 as we 
approach the 'horizon' at r = 0. Note however that this divergence does not effect the volume 
element, ^/—g, and so we should be careful in deciding whether or not these perturbations 
produce a true singularity at r = O.Q In the next section, we will see that in fact with / > 
these solutions are singular, and hence we are only left with a wavy black string for / = in 
which case 

ds'^ = r(f + ^1 du^ + ^dudv + ki{dr^ + r'^dn) (66) 
h \ r J h 

These perturbations represent longitudinal waves carrying momentum along the string with- 
out oscillations. Note that the Qi term in / represents a constant contribution to b{u). 
One may wonder if this wave is really physical or merely an artifact of an awkward choice 
of coordinates, given our theorem on the elusiveness of the wave profile. To see that it is 
indeed physical, we can compute its mass per unit length. Since the oscillating string ([66| ) is 
asymptotically fiat, and we can use the coordinates of eq. (p5D, we can determine the mass 
per unit length according to dE/dy = — (l/Svr) jgEybcdedx^dx'^V^S,^ , where all the tensors are 
defined in the Einstein frame. The vector field C, = dt here is the asymptotic generator of 
time translations, and the integration is carried over a sphere at spatial infinity. The result 
is dE/dy = (Pi + P2 + Q2 + 3(5i + 36(y — t))/6 and since the mass per unit length depends 
on the wave profile, we see that the solution is really a superposition of the string and the 
wave and so is clearly different form the stationary string, where 6 = 0. 

4.1 Parallelly Propagated Singularities 

Having applied the GV technique to the original solution above, we have apparently gen- 
erated a large family of oscillating black string solutions. However, as discussed above, we 
have reason to worry that some of the modes may actually introduce a curvature singularity 
at the null surface which was originally the black string's horizon. Normally the approach 
to proving the existence of a horizon would be to find coordinates in which the metric is 

^The divergence in does indicate a divergence of the norm of the Killing vector h^df^ — 9„, which 
could be interpreted as a geometric singularity |lO) . One should still demonstrate that r = is accessible to 
causal observers, i.e., that this region 'belongs' to the spacetime, as is done in our following analysis. 
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analytic at the null surface in question. For the present undulating solutions, finding such 
coordinates is an enormous problem (see, e.g., 0). While we address this question for the 
monopole waves, i.e., those with / = m = 0, in the following section, the task at hand 
is in fact much simpler. We wish to show that for Z > a given null hypersurface is not 
a horizon, which simply requires finding any geometric quantity which diverges when the 
surface is approached along some geodesic. By the properties of our original black string 
and the approach used to generate their wavy counterparts, the theorem of section ^ tells 
us that no scalar invariants involving the curvature and/or the matter fields contains any 
information about the waves. Therefore given that all scalar invariants are insensitive to the 
oscillations and thus to any singularity which they may introduce, we resort to alternative 
means of probing the wavy geometry. 

Tidal forces prove to be a good tool for resolving our problem. If we approach the null 
surface along a geodesic, and we allow our observer to be slightly non-local {e.g., a string 
which may indeed seem the natural probe in the present context), this observer will be 
able to determine the differences between the gravitational forces acting at different points. 
These forces are determined by the Riemann curvature measured in the rest frame of our 
observer. One may object to the notion of the observer's rest frame, as we have just said 
that the observer of interest is non-local. We will assume that the extension of the observer 
is controllably small, and hence that the center-of-mass frame represents a good reference 
frame in which to express the results. 

Still, identifying a convenient geodesic trajectory to follow in the presence of a general 
oscillation proves to be beyond our abilities. The reason is that one cannot find enough 
integrals of motion to solve the problem in terms of quadratures. We therefore restrict our 
attention to wave profiles which are constant in u. In this case the extra Killing vector, 
h^d^ = du, yields an additional constant of the motion, enabling us to find analytically 
suitable geodesies for any mode. These solutions should be a very good approximation for 
backgrounds with slowly varying wave profiles. When we compute the full curvature of an 
undulating string, we can see that any u dependence only adds contributions of a subleading 
order. Hence if the curvature turns out to be divergent as some hypersurface r = constant is 
approached, the constant profile solutions will contain all the information about the leading 
order of divergences. 

We present our calculations in several steps. First we will show that for each mode of 
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oscillation, there exists a geodesic stretching between the null surface and the asymptotic 
infinity — hence showing that both of these regions belong to the space time. Next we 
will construct the Lorentz transformation relating a natural stationary orthonormal frame 
to the rest frame of the observer moving along the geodesic. Finally we will consider the 
orthonormal frame curvature and boost it to the frame of the infalling observer, in order to 
find the tidal forces which he measures. The divergences found in this way are equivalent 
to parallelly propagated curvature singularities. We will isolate the leading divergences of 
the tidal forces for all modes both at r = and asymptotically, finding that all the localized 
multipoles with / > 1 have unbounded tides on the null surface and that all the growing 
multipole modes with / > 2 have divergences at asymptotic infinity. 

1) Geodesies 

We begin by examining the timelike geodesies of a wavy solution which is excited by 
a single mode. We will demonstrate that there are always geodesies extending between 
r = and the asymptotic region r — oo, and the former is reached in finite affine time 
when starting from finite r. Hence this null surface r = must be included in the manifold 
described by the wavy solution. 

We may write the general solution as 

ds^ = 2F2dudv + F^du^ + (dr"^ + r'^{de'^ + sin^ d d(p^)) (67) 



where 



and ^ = 5(M)r^P^"(cos0)cos(m0 + 5(M)) (68) 



and /, h, k, I are defined in eq. (|5^). To produce a simple real metric, we have expressed the 
angular dependence in terms of an associated Legendre function of the first kind, P[^{cos6), 
as well as the cos{m(f) + 5) factor, rather using spherical harmonics as in eq. (|6lD. Those 
solutions (^) are then reproduced by setting /5 = — (/ + 1) or Z, which corresponds to what 
we will call the localized and the growing modes, respectively. However, for much of the 
following analysis, we will leave this exponent as (3 in order to emphasize the contributions 
coming from the differentiation of this factor. In principle the amplitude B and the phase 5 
are arbitrary functions of m, but as discussed above, to simplify the analysis of the geodesies. 
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we will set both of these to be constants in the following. This will be enough to identify 
the leading divergences, and should still provide good approximation in the case of a slow u 
dependence. In fact we will set 6 = 0, which can be attained with a simple shift of 0. 

As usual to obtain the geodesic equations, we simply consider the Lagrangian L = 
[ds/dX)"^, and write down the Euler-Lagrange equations. Because the Lagrangian does not 
contain a potential, the effective Hamiltonian is conserved, giving {ds/dXy = const. Given 
this integral of motion, we need not consider the Euler-Lagrange equation for the radial 
coordinate. In addition, the two translational Killing vectors, du,dv, produce two more first 
integrals, which leaves us with second order differential equations only for the two angular 
coordinates. Hence the equations for timelike geodesies are 

F2U' = p-uj F2V' + F^u' = p 

2F2v'u' + F^u'^ + ^2 + + sin2 = -1 

^,r^e\, ^r^sin^cos^ ,,2 B ^dPr(cose) , ,2 

2(-^)' = 2 + jT^^^Q ^cosm0 u'^ (69) 

2(^^^^^0')' = -m-^r^Pl^ {cos 9) sinm0 m'^ 

where we expressed the integrals of motion as uo and p. When /3 < 0, one finds that in the 
asymptotic region ( ^5l) these correspond to the energy (which is assumed to be positive), 
and linear momentum along respectively. The last of these equations may be solved by 
setting (f) = (pn = with n = 0,1,..., 2m — 1 for m 7^ 0, while for m = we may 
fix (f) to be any constant. Now with the constant 0, the first term on the RHS of the 
second-to-last equation also vanishes. Hence the latter equation is solved if we now choose 
= 60 = constant where Oq corresponds to an extremum of Legendre function P[^{cos6). 
Then the independent set of equations defining these radial geodesies may be written, using 
F| = (/ + (-)"5r^P;™(cos^o))//i, 



+ {f + {-rBr^Pr{cose,))^ (70) 

'2 F2 

Now we must choose 9q and as well as the constants uj and p, in such a way that our 
geodesies (a) extend to infinity and (b) reach the null surface r = 0. We need the former 



F2 - F2 ' ' ' ' ' ^ hp;' 

r'' ={f + i-TBr^Pricoseo)) - P? + 2^(a; - p)p - F^ (71) 
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to insure that our observer is physically connected to recording devices infinitely far away 
from the gravitational source. The second condition that the geodesic does not turn before 
reaching r = is necessary because we want to probe this region for singularities. 

We begin by considering the localized modes with f3 = —{I + 1) < 0. From the radial 
equation ([fl]), we see that as r oo, for which h, Fi 2 —>■ 1 and — > 0, so r'^ —>■ uj^ —p^ — 1. 
Thus we require u > y/p'^ + 1 so that the geodesies extend to infinity. Now in the limit 
r ^ 0, we have F2 = h^^ oc r, oc and / oc 1/r. Thus the leading contribution on the 
RHS of eq. ( |7TD comes from the term when (3 < —1, and hence we must choose 6q and n 
such that {— )"'BPl^ {cos 6 q) > in order that no turning points occur before reaching r = 0. 
For / > 1 there exist many extrema of P[^{cos6), and it is straightforward to verify that one 
has enough freedom to choose the angles in order that the geodesic reaches r = (even for 
the case m = 0). For the special case (3 = —1 or I = 0, turning points are evaded as long as 
Qi + B> 0.0 

For the growing modes with l3 = I > 0, the analysis is similar to that above. Note that in 
this case, the interpretation of u and p would not be correct since the asymptotic structure 
of the metric is modified. Beginning with the limit r — > 0, we can ignore contribution 
in eq. (^) and in fact one finds that this region is always reached without imposing any 
constraints. In the asymptotic region r 00, the term dominates eq. (|71|) for (3 > 0, 
and we must again choose the angles such that {—)^BP[^ [cos 60) > if the geodesies are 
to extend all the way to infinity without turning back. For the special case j3 = = I, the 
constraint to reach the asymptotic region becomes {uj — p){uj + p + B) > 1 which can always 
be satisfied with an appropriate choice for the integrals of the motion. 

In the following, it will be convenient to absorb (— )"'P;™'(cos^o) into the amplitude. 
Hence we define Bq = (— )"i?P™(cos^o)- For the special case with (3 = —1 and Z = 0, we set 
Bo = Qi + B. 

2) Lorentz Transformation 

Hence we have defined an interesting set of radial geodesies for which the tangent vector 

^Note then that for Z = and B < —Qi, one has solutions for which the radial geodesies can never 
reach r = 0. As pointed out in ref. however, when the y direction is compactified in this case, there 
appear closed timelike curves, which can communicate with observers at infinity, implying the breakdown of 
chronology. In any event, we will ignore this pathological case in the following. 
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is given by = dx^ /d\ where 




r' = ± ( (/ + B„r^)-i^{uj - pf + 2^(i^ - p)p - F'i^ ' (72) 

and 9' = Q = cf)' . Here the — (+) sign corresponds to inward (outward) directed geodesies. 
We will want to examine the tidal forces in the rest frame of an observer moving with this 
five- velocity. Hence as an intermediate step, we determine the Lorentz transformation which 
takes us from a natural stationary frame, in which the curvature is easily calculated, to the 
observer's freely-falling frame. 

First to define our stationary orthonormal frame, we complete the squares in our general 
metric ( |67D 

ds^ = -^dv^ + Flidu" + ji^dvf + (dr^ + r\d9^ + sin^ 9 dcj)^)) (73) 

We see an obvious orthonormal basis of one-forms is 

^2, 4 r , , ^2, r dr g rd9 ^ rsm9d(f) 

e = —dv e = t^du + —dv e = — e = —— = (74) 

F3 F3 Fi Fi Fi 

Note that in this basis, e° is distinguished as the unit time-like one-form, at least everywhere 
along our radial geodesics.[| In this frame, we have = e°'^V^ where the e"^ are the 



components of the fiinfbein ([7^) 



^3. . , P P T/r , f ^3 r n2 , ' P)P 



v" = yi--p) + y y' = Y v^^ = ± t4(--p)^ + 2"^-i (75) 

^2 ^3 ^3 \^2 ^2 ) 

along with = = . As a check, one may easily verify that rjabV^V^ = ~1- 

Now we wish to find a Lorentz transformation which takes a unit time-like vector A^" = 6q 

into the observer's five-velocity: V°- = L^-^N^. Then applying this transformation to our 

stationary fiinfbein ( |71D would produce a natural basis of orthonormal one-forms which the 

observer might use in his rest frame. Of course, we are left with some ambiguity in the choice 

^This is guaranteed since = is trivially positive, while we ensured that F| > in order that the 
geodesies reached between r = and r ^ oo. The sole exception, which we ignore, is the case (3 = = 1 for 
which we could have < if i? < — 1. 
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of the S0(l,4) matrix defining this Lorentz transformation. One approach to resolving this 
ambiguity is defining the transformation by parallelly propagating the stationary frame in 
from infinity along our radial geodesic. Since = = V^, a much less labor-intensive 
approach is to simply take the uniquely defined SO (1,2) matrix mixing only the 0, 4 and r 
directions. The latter will differ from that defined through parallel propagation by an SO (4) 
transformation mixing the spatial one-forms. Such a rotation however will not introduce any 
new divergences, and hence it suffices to consider the simpler boost. One can think of this 



transformation as the result of parallelly propagating in not our stationary fiinfbein ( fM]) but 
rather a rotated version of it. The simpler S0(l,2) Lorentz transformation may be written 

as 



yr 



V 



1 + 



yO(v'0+i) 

yi yr 



yO(yO + l) 






1 + 



yr 

y4 yr 

yo(yo+i) 



VO{vo+i) 











1 










1/ 



(76) 



At this point let us comment at the behavior of the five- velocity in the regions of interest. 
For P < —1 and r — > 0, the dominant factor is oc r^^^^^^^"^ leading to ^ (cu — 
p)\^BqQ2 r^'^^^)/^ ~ V while ~ 0. Hence the observer is accelerated to almost a null 
radial geodesic as he nears r = 0, and L approaches an infinite boost in the radial direction. 
For /? > 0, F3 and all of the components of are finite at this null surface. Similarly 
for P > and r 00, F3 oc r^/^ again dominates to producing the five-velocity ~ 



(cu — p)^/BQr^^'^ ~ and V^^ ~ 0. Hence eq. (|76|) yields an infinite radial boost as r 00. 
Finally, we note that for 13 <1 and /3 = 0, F3 and all of the components of V"' and L remain 
finite in the asymptotic region. 

3) Divergent Tides 

Now we wish to examine the tidal forces experienced by an observer following our radial 
geodesies. We will focus here on the curvature components with all indices in the 0, 4, r 
subspace, and study the corresponding components in the observer's rest frame. Since our 
Lorentz transformation (^) does not mix these indices with the 9, cj) directions, it simpli- 
fies our presentation to only consider these components. Our omission of the components 
carrying angle indices does not imply that they are finite. In fact the components with 
two indices in the 9, cj) subspace also typically diverge whenever we find divergences in the 
following analysis, and in qualitatively the same way, as discussed at the end of this section. 
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Our final results are that divergent gravitational tides appear on the null surface r = for 
the localized modes with (3 < —1 and at asymptotic infinity for the growing modes with 
f3>2. 

With our definition of the orthonormal basis ([74D, it is straightforward to compute the 
curvature in this frame. As discussed above, we focus on the curvatures in the 0, 4, r subspace 
for which the nonzero components are 

D4r4r _ p' F' I ^ F'F' ^P" ^ ( F')"^ ^ ( F')'^ 

^ - f/'^'^ F^f/'^' f/' Ft^'' 4Ft^'> 



P p2 p2 p2 

pOrOr _ _f2 p' p' _L 1 P' _ 11 _ llf P'^2 

T 2 ~r T-1 ^ 

2 -^2 



+ 11 + lliT' f: - (77) 



T-,0404 jIL /p/^2 



= -^F;,F[ + -^F;,F^-fF;;-^^{F;,r 

^3 ^2^3 ^3 ^3 



p2 p p2 

.lJ_p" + —F'F' - —^(F')^ 



where the primes indicate partial derivatives with respect to r. 

Now we wish to consider r 0. To simplify our calculations further, we only consider 
the leading order of these curvature components. This approach is consistent, because we 
will see that all of the nonzero components ( [77D are of the same order of magnitude at the 
null surface. In fact these terms are all finite there, and so any divergence in the tidal 
forces can only arise from the boost matrix (|76| ) upon transforming to the observer's rest 
frame. Given the discussion at the end of the previous subsection then, it is clear that the 
only possible divergences will occur for the modes (3 = —{I + 1) < 0. For these modes, 
as r ^ 0, Fi ^ r/^/lW2, F[ ^ 1/v^, ^ -(Pi + P^) / [P.P^y^'^ F^ ^ r/Q^, 
F^ ^ I/Q2, F^' ^ -2/Q2 and F, ^ jBo/Q,r-'/\ F^ ^ ~{l l2)JB^^r-''^~' and F^' ^ 



(/(/ + 2,) /A)^ Bq/ ^, and so the above curvature components (fT^) reduce to 
^4,4, ^ l + 2/j(/3+l) ^0,0, i_2/?(/? + l) 



4P1P2 4P1P2 



d0404 ^ D4r0r ^ /^(/^ + 1) t^Q\ 
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Transforming to the observer's rest frame, R"-^'^'^ = L"" kV' iL'^ ^L'^ ^R^^^^ ■, we find 

1 

4P1P2 



^ahcd _^ _J_^^abcd _ 2/3(/3 + i)Af "'^j (79) 



where 



^ 1 Ar^ = o (80) 



^0404 ^ {y^f^[u;-pfB,Q,r^-^ At'^ ^V\vr^p{u;-pf^BoQlr^^-'^/' 

Hence we see that there are divergences but that they only appear in the terms proportional 
to /?(/5 + 1). The /^-independent contributions are essentially boost invariant. These diver- 
gences are therefore absent for /? = — 1 or / = 0. This result might have been expected since 
in this case, for a constant wave profile, the solution is essentially the original regular black 
string with a modified Qi. The leading divergences are easily seen to be 

A4.4r ^ l{l + l)BoQ,{u-pf ^0,0, ^ 1 - 2/(Z + 1) 

2 PiP2r^+^ 4P1P2 



A0404 Ijl + 1) BoQ^jio - pf £,4,0, l{l + 1) ^ BqQUu^ ~ pfp 

^ 2 PiP2r'+2 2 P^P^r^/2+^ 

So we see that for all higher multipoles with (3 = —(/+!) < — 1 or / > 0, there appear singular 
tidal forces on the null surface r = 0. Because these divergences will not be cancelled by any 
other terms of the metric for slowly oscillating strings, we conclude that all these space-times 
have a null singularity at r = 0. Hence, the excitation of these higher modes on the black 
string results in the appearance of naked singularities, and so the resulting solutions are no 
longer black strings after all. Thus we rule out all of these higher multipoles as a variety of 
non-stationary hair. 

A similar calculation shows that the growing wave modes with (3 = I > 1 have diverging 
tidal forces in the limit r oo. In this limit. Pi 1, F[, F" — > 0, P2 — 1, P2, P2' and 

P3 ^ v^r^/^ P;^ ^ (/3/2)/B;;r^/2-i and P^' ^ {(3{(3 - 2)/4)v^r^/2-2^ ^lence to the 



leading order, the curvature components (^) become 



2^2 
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The boosted curvatures, to the leading order, are given by 

The terms proportional to A"^'^'^ are all of the subleading order. The limiting values of /S.'^'^'^ 
when r — > cxD are given by 

^f^r _^ _ pfB^rl^ A°'^°'- ^ 1 (84) 

and so the leading divergences are easily seen to be 

^4r.0r _^ - p + p(^ - p)2)/B;r^/2-2 (85) 

Hence we find diverging tides in the asymptotic region for all /5 > 2. For /3 = 2 we find 
finite tides at infinity; however, this implies that at infinity the energy density approaches a 
constant, and hence the total energy per unit length of this wave diverges. Indeed, this might 
have been expected as these solutions are not asymptotically flat. Instead, they represent 
geometries with gravitational wave energy concentrated far away from the black string. 

To close this section, we will discuss the results for the remaining orthonormal components 
of the curvature tensor. It turns out that there are two cases to consider: R'^^'^'^ and R"-'^^!^ 
where a, (3 take values in 0, 4, r while the remaining frame indices are 9 or 0. Straightforward 
evaluation shows that the curvature components with an odd number of angle indices vanish 
when evaluated on our radial geodesies which were chosen so that dgF^ = = d^F^. Further, 
as the components in eq. (ffTl), the nonvanishing components considered here remain finite in 
the stationary frame (|7iD. So clearly the boost (|76[), which is trivial in the angle directions, 
cannot introduce any divergences in R"'^'^'^. Now for the i?"^"^^, we have symbolically R = L^R 
upon boosting to the observer's rest frame, where we have only indicated the non-trivial 
components of the boost matrix. Because the are everywhere finite, the tidal forces 

here diverge only as badly as the leading order divergence in L^. As r — 0, this is l/r'+^ 
for P = —(/ + !), and as r ^ oo, this is for [3 = 1, just as found above. Note that 
various cancellations above reduced the singularities from what one might expect for L^. 
Also as above, the original curvature components supply a factor of /(/ + 1) multiplying 
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these divergent terms so that the special cases j3 = —1,0,1 survive without any divergent 
tidal forces. These qualitative arguments are confirmed with direct evaluation. Hence one 
finds that there are no worse divergences than found by considering only the 0, 4, r subspace. 

Again, our final conclusion is that all of the modes with /3 = — (Z + l) < — lin fact produce 
null singularities at r = 0, while for (3 = I > 1, singularities appear in the asymptotic region 
r — s> oo. 



5 Longitudinal and Transverse Oscillations 



We now return to consider the case of the longitudinal waves ( |5BD with / = 0, = — 1 and that 
of the transverse waves (|6^) with I = 1 = (3. The preceding analysis revealed no divergent 
tidal forces in either of these cases. There is still the possibility that divergent tidals might 
occur when the analysis is extended to consider derivatives of the curvature (much like the 



case discussed in [|23[). In the case of the transverse waves, the worry would be that problems 



arise in the asymptotic region where the Lorentz boost (|7q) is divergent. However, recall 



that a coordinate transformation (|6^) was found for which the metric became manifestly 

asymptotically flat.0 Hence one will never find any divergent tidal forces as r — > oo. In fact, 

for the transversal waves one can argue that the divergence of the Lorentz transformation 

([f6|) is not physical, but comes from an incorrect choice of gauge in the limit r oo. Namely, 

in this case one could view the asymptotically flat form of the transverse wave metric ( |65D 

with Ai = const, as the correct example of the leading order behavior of slowly oscillating 

transverse waves. The fact that this solution is asymptotically fiat implies that the Lorentz 

transformation ([76|) must be finite in the limit r — oo. 

For the longitudinal waves, the potential problem would be at r = where again divergent 

components appear in the transformation ([7B|). A definitive demonstration of the regularity 

of these solutions would require finding a coordinate transformation for which the metric ( pBD 

becomes analytic at the null surface. For a wave profile constant in u, we recover the original 

solution (^) with a modified charge in which case it is straightforward to find the analogue 

of Eddington-Finkelstein coordinates, both past and future. For the nonstationary case, 

Horowitz and Marolf |§ have shown that this solution has a continuous (but not necessarily 

^Recall that no divergent tidals were found for ^ = = /? either, but in eq. (^) this mode was shown to 
be purely the result of a coordinate transformation. 
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smooth) metric at r = 0. Further, our theorem shows that all the scalar curvature invariants 
of this solution are identical to the original black string, while our preceding analysis found 
no evidence of diverging tidal forces for slowly oscillating strings. In fact, we can make 
an even stronger statement about the 'invisibility' of the monopole wave: if we compute 
the curvature of the oscillating string in an appropriate orthonormal basis (see below), we 
find that it does not depend on the wave at all! Given all these results, one is tempted to 
conjecture that the longitudinal waves are completely regular, and so we are lead to attempt 
a construction of analytic coordinates. It turns out that we can find such coordinates at 
least for a certain class of wave profiles. A rather surprising result in our analysis is the 
quantization of a certain constant present in the test function, coming from the requirement 
that all the derivatives of the metric be continuous as the null surface is approached. It is 
not clear to us at present whether this quantization is an artifact of our ansatz or whether 
it really represents a physical effect. We will defer a more detailed investigation of this issue 
to future work. At this point, we can only provide some guidelines explaining how we have 
arrived at the analytic ansatz and the associated quantization condition. 

We begin by considering the longitudinal waves described by the metric (|66D which may 
be written: 



ds'={l + ^' + ^^""A du' + -^dudv + (r + P{){t + P,){'^ + dn\ (86) 
- ' Q2 r + Q2 \ H ' 



Now the coordinate transformation, dv = dv + ^^—^^^^^du, simplifies the metric somewhat 
producing^ 

2r f dr'^ \ 

ds' = p\u) du' + —dudv + (r + P{){r + P2) [ ^ + d^l] (87) 

r + Q2 I 



where 



2/ N Qi + b{u) 
P W = 



Q2 

Implicitly we assume that p'{u) > — see footnote 6 — and as usual dQ is the metric on a 

two-sphere. We see that the location of the null surface remains at r = and that it appears 

to be a metric singularity: we know however, that if p' were a constant, we could easily 

show that this singularity is just a coordinate artifact, as we have mentioned above. Let us 

^If one now constructs an orthonormal frame analogous to that ( [t^ found for the metric (^3|), one would 
find that the wave profile p^{u) completely disappears from the frame components of the curvature tensor. 
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therefore try to follow this argument as closely as possible. We can put the {v, u) part of 
the metric in the form conformal to the constant case, by defining a new coordinate z: 
dz = p^{u)du and q{z) = l/p{u) in which case the metric becomes 



ds' 



q\z) 



dz^ + 



2r 



r + Q: 



-dz dv 



dr"^ 

+ {r + Pi){r + P2) \ — + dn 



(89) 



Now we can introduce new coordinates v and z in analogy with the tortoise coordinates we 
would usually define for a stationary solution. A bit of algebra leads to the choice 



dv 



dz 



dv + 



Pi + P2dr, 2P1P2 
r + Q2)[r + 



P1 + P2' 



[Podr , „ , 
V Pi r 



(90) 



These coordinates are oriented towards the future portion of the null surface, and they 
clearly show the location of the future horizon as r ^ along with v +00 and z —00. 
Reversing the signs of the shifts to v and z, we can go to the past horizon. These coordinate 
transformations are designed to 'eat up' the metric singularity manifest in the factor 
appearing in Qrr- If we rewrite the metric (|88|) in terms of the tortoise coordinates (|90D, we 
find 



ds' 



q\z[ 



dz^ + 2- 



L r + Q2 

+ (r + Pi)(r + P2 



dz dv + 2\l -TT- — ^ dr dv 



Pir + Q2 



P1-P2 
VPJ^ 



dr dz 



1 - q\z) 



dr'^ + dVL 



(91) 



Here z is to be understood as an implicit function of z and r: z = z + + -Pi log('"/-Pi)], 

and it diverges to —00 as r ^ as dictated by the logarithm. Now if = 1, the divergent Qrr 
term in (0) would be absent, and we would obtain a smooth metric at the null surface r = 0. 
Can we now select a wave profile q^{z) such that a similar cancellation still occurs in the limit 
r — s> 0? The answer turns out to be in the affirmative: Consider a function q^{z) which in the 
limit z —00 converges to 1 — A^exp((2 + a)z / ^ P1P2) for some positive a. Substituting 
the appropriate coordinate transformation z = z{z,r) in this expression, we get that in the 
limit r the wave profile to the leading order is = 1 — A^(r/Pi)^+"e'^*^^'''\ where for 
convenience we have defined the linear function F{z, r) = (2 + a){z + \J P2/ Pi r) / ^/PiPi- As 
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long as a > 0, this is precisely of the form needed to cancel the pole in QrJ The metric, in 
this limit, becomes 



which is smooth at r = 0. 

At this point, however, we can push these arguments one step further. A glance at the 



coordinates z and t are well defined as the horizon is approached — we ignore subleading 
terms in for the moment. Moreover, if a is chosen to be an non-negative integer, i.e., 
a = 0, 1, 2, ... , it is easily verified that all the derivatives with respect to r are also well 
defined in this limit. Indeed, we see that the only possibly contentious terms in the metric 
are the factors r^"*"" and r". If a were not integral with > a > — 1, then taking A^ 
r-derivatives of the metric would produce factors which diverge as r 0. At this point, we 
cannot tell whether this divergence could cause some covariant derivative of the curvature 
tensor to diverge as well (resulting in a null singularity as in ||2^), or if it could be removed 
by another even more clever change of coordinates. On the other hand for integral a, we see 
that after taking some number of derivatives of the metric with respect to r, the contentious 
factors disappear altogether, leaving an expression which is perfectly well defined as r — 
and z — > — oo. 

It is a matter of a simple counting of powers to convince oneself that all these conclusions 
remain unchanged if the function q'^{z) can be written as a uniformly convergent power 
series q'^{z) = 1 — Z]^2^n'"^" = exp{nz / P1P2) and z —00. The subleading 

contributions would be integer powers of the leading term, therefore still having well-defined 
r-derivatives. The requirement of uniform convergence ensures that the summations and 
derivatives commute, resulting in the conclusion that the profile is itself analytic on 
the horizon. This allows us to extend our arguments in order to analytically continue the 
solution through the past horizon z — > 00, as well as the future horizon at z ^ —00. For 
example, consider a wave profile of the form q'^{z) = 1 — A^/ cosh.{kz / \^ P1P2) for some 




(92) 



metric ( p2D immediately shows that all of its derivatives with respect to the angles and the 
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positive integer > 2.|^ Expanding this profile in the hmit z 



oo (future horizon), we 



find g^(z) = \ — X]m=o(~l)"' 6^P((2"^ + l)A;z/-\/PiP2), precisely of the form guaranteeing 
the existence of the limit r ^ 0. If we approach the past horizon instead, with z —>■ oo, 
we find q'^{z) = 1 — Z]m=o(~l)™ ^^P("~(2"^ + 1)^-2/^-^1-^2), again rendering the limit 
well-defined. Hence we see that with this example we really have constructed a wavy black 
string with regular future and past horizons. Of course, this example is easily generalized 
e.g., by taking linear combinations of inverse cosh's with different integers k, or by taking 
a wave profile of the form A'^/(w'^ + w''^') with independent k and k'. Hence as we have 
claimed above, we see that at least for a certain class of wave profiles, the longitudinal wave 
solutions are analytic at r = 0, which therefore can be identified as a regular event horizon. 
Hence we can think of these waves as time-dependent hair on the black string. Clearly, our 
argument is only an existence proof. It would be interesting to prove that the general family 
of longitudinal waves is regular, or to determine the precise conditions which the wave profile 
must satisfy in order to ensure regularity. We note here that for all of our examples above the 
profiles have essentially compact support in the z coordinate, and therefore in the horizon 
limit the waves are exponentially damped, the solution approaching the stationary string. 

One could also consider introducing transverse and longitudinal waves on the black string 
at the same time. This would amount to adding to the line element (|86|) a term of the form. 



In this case, we would perform all of the same transformations described above with the 



where we have rewritten a{u) as a{z). Now regularity of the horizon is secured by demanding 
that a{z) is an analytic function of z and r. The latter is easily accomplished with wave 
forms similar to those discussed for q{z)'^ above. 



Note that we would require < 1 in order to ensure that = > 0, as assumed above. 



e.g. 




(93) 



result that the final metric (|9TD is modified by the addition of 




(94) 
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6 Discussion 



In this work, we have investigated geometric properties of Garfinkle-Vachaspati waves, ob- 
tained by superposing gravitational oscillations on stationary solutions with null hypersurface- 
orthogonal Killing vectors. We have first given a detailed discussion of the adaptation of 
this solution generating technique to a variety of supergravity models, relaxing some of the 
matter sector constraints imposed previously |T^. Then we have developed a purely geo- 
metric theorem, stating that the GV modes are completely invisible to any scalar invariants 
constructed from the metric and matter. The proof of our theorem does not rely on the 
dynamics at all. It is a consequence of the null symmetry, which must be present in order to 
interpret the solutions as waves, and hence it applies to any metric which can be represented 
in the generalized Kerr-Schild form with respect to the symmetry. Thus our theorem holds 
in any theory of gravity and in an arbitrary number of dimensions. 

As a consequence, it is evident that one must further scrutinize wavy string geometries 
using non-invariant probes, such as tidal forces. Recall that we have in fact shown that none 
of the scalar invariants changes under the Garfinkle-Vachaspati map. In the group-theoretic 
language, they are invariants of the orbits of the solution-generating map. Now, this could 
suggest that the wavy strings might contain curvature singularities identical to those of the 
corresponding stationary solutions belonging to the same orbit. Consider, e.g., the stationary 
solutions with all charges equal. Since they correspond to the extremal Reissner-Nordstrom 
black strings, with traversible horizons and spacelike curvature singularities at the core, one 
could be tempted to conjecture that all the GV wavy modes superposed on such solutions 
also have similar curvature singularities - for the curvature scalars are exactly the same, and 
in particular also blow up at the 'core'. However, in general this is incorrect. The argument 
is incomplete because it does not answer whether such singular regions can be causally 
connected to the physically interesting parts of the spacetime, such as the asymptotic infinity. 
Thus it becomes evident that the notion of singularity is a far more complex issue here: in 
the examples of oscillating solutions we have encountered, some of the singularities encoded 
in scalar invariants may not belong to the physical spacetime, because no causal geodesies 
can ever reach it. The simplest example of such behavior is given by the solutions discussed 
by Gibbons, Horowitz and Townsend p2[. These solutions correspond to the choice of the 
wave degree of freedom such that / -|- h{u)/r = 1 in eq. (^). Although they belong 
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to the same GV orbit as the stationary extremal Reissner-Nordstrom strings, their causal 
structure is dramatically different: in order to get the analytic extension across the horizon, 
one has to extricate the region normally associated with the 'Interior' of the black string from 
the physical spacetime, and replace it with a mirror image of the 'Exterior'. The resulting 
manifold is a completely nonsingular spacetime with a horizon, with the topology of the 
Anti-de-Sitter horizon crossed with a sphere, dividing two asymptotically fiat regions. Non- 
curvature singularities, manifest in the presence of geodesic incompleteness, may reemerge 
if one attempts to compactify some of the longitudinal directions in order to descend to 
lower dimensions. Nevertheless these singularities are typically milder than the original 
ones which got cut out in the process of analytic extension. Therefore we can see that our 
theorem states that scalar invariants are identical as analytic functions of the coordinates 
on each GV orbit. The observers 'surfing' on GV waves of differing geometries however may 
choose to access different regions of this 'complex' coordinate plane as physical, and hence 
see different divergences. We believe that this unusual behavior is peculiar to wavy strings 
(or more generally wavy p-branes) and hence pathological. 

With this in hindsight, we have found that most of the undulating solutions cannot be 
interpreted as black strings: the excitation of any localized multipole mode with / > 1 gives 
rise to a naked null singularity as opposed to an event horizon. Similarly, oscillations of 
growing modes with / > 2 result in divergences very far away from the string. To uncover 
these divergences, we have found a timelike geodesic for each excited multipole, which extends 
to infinity and reaches the null surface in a finite affine time. We have then propagated an 
observer along this geodesic to the vicinity of the region of geometry we wished to explore, 
and isolated the leading divergences of the curvature forms Lorentz-transformed to the rest- 
frame of the infalling observer. These divergences translate into infinite physical forces 
detectable by the observer, and hence must be interpreted as naked singularities. 

One could easily generalize these conclusions to include the waves propagating in the 
remaining 'internal' dimensions, which we have ignored until now because they were passive 
spectators for the most part of our investigation. In the framework of heterotic string theory 
this can be done as follows: we could take the tensor product of the string frame solution 
(^) with a four-torus, lifting it to nine dimensions. This is guaranteed to be a solution 
of the effective action of heterotic string theory in nine dimensions, with the matter fields 



identical to those of the original five-dimensional solution (|52D. We could then view the 
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Laplacian constraint of (|T8D as a linear combination of the spacetime and internal parts. 
From the solution (^2]), the Laplacian constraint in the Einstein frame in nine dimensions 
can be decomposed as V^\l/ + /c/V^\E' = 0, with k and / defined in eq. (^41). The operators 
and Vy denote the flat space Laplacians in the spatial sections and the internal 
space, respectively. Because the internal space is a four-torus, we can expand the wave 
profile \& in the Fourier series with respect to the internal space basis, which consists of the 
exponentials of linear combinations of the internal coordinates exp(c ■ y). The coefficients 
of the coordinates are integer multiples of the inverse radii of the torus, and would give rise 
to the mass of the transversal part of the wave profile = J2i=iC^- In particular, the 
zero modes, corresponding to Ci = ■ ■ ■ = C4 = would trivially reduce to the case of the 
massless spacetime waves we have discussed in the main part of the paper. These are the 
only scalar harmonics on the torus and hence we can ignore them. All the massive modes 
however turn out to be naked singularities. This can be seen as follows: the spacetime part 
of the Laplacian for each massive mode of the wave profile becomes V^\I/m = M'^kl'^M- 
Now using the spherical symmetry of the stationary solution to implement the separation of 
variables = RMiYirm we at last obtain the equation for the radial degree of freedom: 

t^Rmi + 2rRMi - + 1) + M^PiPa + M\P^ + P2)r + M\^)Rmi = (95) 

From the general theory of differential equations, we know that this equation has two classes 
of solutions - localized (approaching a constant as r ^ 00) and growing (diverging in the 
limit r 00). The latter do not represent wavy strings but rather waves filling up the whole 
spacetime — and hence lead to divergences far away from the string, as in the case of the 
massless growing modes discussed in section 4. This can be seen as follows: substituting 
Rmi = T^Mi/f in ( p5D and keeping only the leading terms in the limit r 00, we find 
TZmi = M'^TZmi- Hence the growing modes diverge at infinity as Rmi — exp{Mr)/r — 
i.e., faster than any power, and so they must produce naked singularities far away from 
the string for all /. The former case is more interesting, because by the same arguments as 
above, far away from the string these solutions damp out as exp(— Mr)/r. Thus they might 
be interpreted as wavy strings — except that they are singular in the limit r — > 0. To show 
this, we note that the point r = is a regular singular point of the differential equation (|95|) . 
Hence all the solutions admit Laurent series representation, with the leading power being r^. 
From dH) we find /3 = (-1 ± ^1 + Al{l + 1) + AM^PiPi) /2. The plus sign in the definition 
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of /3 leads to finite (in fact, vanishing) limits of the radial function as r ^ and, as we know 
from the theory of differential equations, it must correspond to the growing solutions which 
are singular at infinity. The minus sign gives solutions which diverge as Rmi — l/r'^' on 
the null surface. In particular, even when / = 0, the wave profile diverges with the negative 



power j3 = — (1 + Vl + 4:M'^PiP2) /2. As we have seen in section 4, such divergences could 
produce infinite tides — and since now + 1) = M^PiP2 7^ 0, the tides we have computed 
there diverge even for the monopole mode! Similarly, we can see that other localized massive 
multipole modes also have divergent tides on the null surface. Hence we see that all the 
nontrivial internal modes lead to naked singularities, which can be detected by tidal forces. 

In the case of localized / = and growing / = 1 modes, we have found a special subclass of 
wavy strings which are analytic on, and everywhere outside of, the null surface, and not only 
continuous - thus extending the earlier results of Horowitz and Marolf |TD[. A puzzling 



feature of our special family of analytic black strings with non- stationary hair is obviously 
the quantized nature of the parameter a. The quantization of a in a technical sense is more 
similar to, say, the quantization of the radial modes of an electron in a hydrogen atom, 
where we also require that the wave function of the system is analytic — yet the obvious 
absence of any quantum-mechanical scales precludes the possibility that our quantization is 
of microscopic origin. On the other hand, it is possible that this is just an artifact of our 
choice of the form of the wave profile yielding analytic metric. In our opinion, this ambiguity 
merits further interest. 

Finally, we should remark once more that our analytic examples comprise only an exis- 
tence proof of time- dependent hair. They do not encompass the general / = and / = 1 
modes. Specifically, we do not cover arguably the most interesting case when the coordinate 
along the string is compact, which is a candidate for a hairy black hole in four dimensions. 
Given our demonstration above that the excitation of the passive internal dimensions leads 
to naked singularities, we should note here that the compact wavy strings may also be sin- 
gular. A simple argument supporting this could be based on the fact that the compact wave 
profiles are periodic functions. Hence the exponential damping of the wave in the horizon 
limit, manifest in eq. (0) and crucial for demonstrating analyticity of the metric in our 
examples, is absent. This actually turns out to be correct — but to prove that the compact 
wavy strings are indeed singular requires a careful examination of the tides beyond the lead- 
ing order, and will be given in the forthcoming work of Horowitz and Yang pSl . So, while 
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our examples show that it is possible to have wavy strings with regular horizons, and that 
the waves can be understood as time-dependent hair on black objects, this identification 
must be applied sparingly and discriminatively. In general, the wavy solutions are not hairy 
black strings, or decompactified black holes, because they contain regions with infinite tides, 
which render the solutions singular. Thus it appears unlikely that the classical wave modes 
alone can give the complete accounting of the microscopic black hole degrees of freedom, as 
proposed in |Tl|, |I5|. The excitation of most of these modes turns the original black string 
into a very bright object. Hence such modes cannot be interpreted as black hole hair - in a 
manner of speaking, they are too curly to weave a smooth spacetime. 
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